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Abstract 

The long-range Frohlich electron-phonon interaction has been identified as the most essential 
for pairing in high-temperature superconductors owing to poor screening, as is now confirmed by 
optical, isotope substitution, recent photoemission and some other measurements. I argue that low 
energy physics in cuprate superconductors is that of superlight small bipolarons, which are real- 
space hole pairs dressed by phonons in doped charge-transfer Mott insulators. They are itinerant 
quasiparticles existing in the Bloch states at low temperatures as also confirmed by continuous-time 
quantum Monte-Carlo algorithm (CTQMC) fully taking into account realistic Coulomb and long- 
range Frohlich interactions. Here I suggest that a parameter-free evaluation of T c , unusual upper 
critical fields, the normal state Nernst effect, diamagnetism, the Hall-Lorenz numbers and giant 
proximity effects strongly support the three-dimensional (3D) Bose-Einstein condensation of mobile 
small bipolarons with zero off-diagonal order parameter above the resistive critical temperature T c 
at variance with phase fluctuation scenarios of cuprates. 

PACS numbers: 71.38.-k, 74.40. +k, 72.15.Jf, 74.72.-h, 74.25.Fy 



I. INTRODUCTION: ESSENTIAL PAIRING INTERACTION IN CUPRATES 



Although high-temperature superconductivity (HTS) has not yet been targeted as 'the 
shame and despair of theoretical physics', - a label attributed to low-temperature supercon- 
ductivity during the first half-century after its discovery - controversy of current theoretical 
constructions has led many researchers to say that there is no theory of HTS and no progress 
in understanding the phenomenon. A significant fraction of theoretical research in the field 
has suggested that the interaction in novel superconductors is essentially repulsive and un- 
retarted, and it could provide high T c without phonons. Indeed strong onsite repulsive 
correlations (Hubbard U) are essential in shaping the insulating state of undoped (parent) 
compounds. Different from conventional band-structure insulators with completely filled 
and empty Bloch bands, the Mott insulator arises from a potentially metallic half-filled 
band as a result of the Coulomb blockade of electron tunnelling to neighboring sites l|. 
However, the Hubbard U model shares an inherent difficulty in determining the order when 
the Mott-Hubbard insulator is doped. While some groups have claimed that it describes 
high-T c superconductivity at finite doping, other authors could not find any superconducting 
instability. Therefore it has been concluded that models of this kind are highly conflicting 
and confuse the issue by exaggerating the magnetism rather than clarifying it [2J. 

Here I discuss a multi-polaron approach to the problem based on the bipolaronic exten- 
sion of the BCS theory to the strong-coupling regime |3j. Attractive electron correlations, 
prerequisite to any superconductivity, are caused by an almost unretarted electron-phonon 
(e-ph) interaction sufficient to overcome the direct Coulomb repulsion in this regime. Low 
energy physics is that of small polarons and bipolarons, which are real-space electron (hole) 
pairs dressed by phonons. They are itinerant quasiparticles existing in the Bloch states at 
temperatures below the characteristic phonon frequency. Since there is almost no retarda- 
tion (i.e. no Tolmachev-Morel-Anderson logarithm) reducing the Coulomb repulsion, e-ph 
interactions should be relatively strong to overcome the direct Coulomb repulsion, so carriers 
must be polaronic to form pairs in novel superconductors. 

In our approach to cuprate superconductors we take the view that cuprates and related 
transition metal oxides are charge-transfer Mott-Hubbard insulators at any relevant level of 

n 

doping [3(. The one-particle density-of-states (DOS) of cuprates is schematically represented 
by Fig.l, as it has been established in a number of site- selective experiments [4J] and in 
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FIG. 1: DOS in cuprates. The chemical potential \x is inside the charge transfer gap as observed 



in the tunnelling experiments 



3] because of bipolaron formation [3] . It could enter the oxygen 



band in overdoped cuprates, if bipolarons coexist with unpaired degenerate polarons 
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the first-principle numerical ("LDA+U") 5] and semi-analytical cluster [6J band structure 
calculations properly taking into account the strong on-site repulsion. 

Here d-band of the transition metal (Cu) is split into the lower and upper Hubbard 
bands by the on-site repulsive interaction U, while the first band to be doped is an oxygen 
band within the Hubbard gap. The oxygen band is completely filled in parent insulators, 
and doped p-holes interact with phonons and with spin fluctuations of d-band electrons. 
A characteristic magnetic interaction, which could be responsible for pairing is the spin- 
exchange interaction, J = 4t 2 /U, of the order of 0.1 eV (here t is the hopping integral). On 
the other hand, a simple parameter-free estimate of the Frohlich electron-phonon interaction 
(routinely neglected within the Hubbard U and/or t — J approach) yields the effective 
attraction as high as 1 eV {jj. This estimate is obtained using the familiar expression for 
the po.axon .eve. shift, E r , the hi g h-f_ ^, and the static, e„, dielectric constants of 
the host insulator, measured experimentally [9j, 



where n^ 1 = e^ 1 — Cq 1 and the size of the integration region is the Brillouin zone (BZ). Since 
€oo = 5 and eo = 30 in L^CuC^ one obtains E p = 0.65 eV. Hence the attraction, which is 
about 2E p , induced by the long-range lattice deformation in parent cuprates is one order of 
magnitude larger than the exchange magnetic interaction. There is virtually no screening of 
e-ph interactions with c— axis polarized optical phonons in doped cuprates because the upper 
limit for an out-of-plane plasmon frequency (< 200 cm -1 )[l0j is well below characteristic 
phonon frequencies, u ~ 400 - 1000 cm _1 . Hence the Frohlich interaction remains the most 
essential pairing interaction at any doping. 

Further compelling evidence for the strong e-ph interaction has come from isotope ef- 
ects [ill], more recent high resolution angle resolved photo emission spectroscopies (ARPES) 
lgj , and a number of earlier optical [13|, |l4j, |l5|, |l6j and neutron-scattering [l7| studies 



of cuprates. The strong coupling with optical phonons, unambiguously established in all 
high-temperature superconductors, transforms holes into lattice mobile polarons and mobile 
superconducting bipolarons as has been proposed 
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When the e-ph interaction binds holes into intersite oxygen bipolarons [3j, the chemical 
potential remains pinned inside the charge transfer gap. It is found at a half of the bipolaron 
binding energy, Fig.l, above the oxygen band edge shifted by the polaron level shift E p , 
as clearly observed in the tunnelling experiments by Bozovic et al. in optimally doped 
Lai. 85 Sr 15 Cu0 4 ]7]. The bipolaron binding energy as well as the singlet-triplet bipolaron 
exchange energy (section 3) are thought to be the orig in of normal state charge and spin 



pseudogaps, respectively, as has been proposed by us [22j and later found experimentally [23 ]. 
In overdoped samples carriers screen part of the e-ph interaction with low frequency phonons. 
Hence, the bipolaron binding energy decreases 24| and the hole bandwidth increases with 
doping. As a result, the chemical potential could enter the oxygen band in overdoped samples 
because of an overlap of the bipolaron and polaron bands, so a Fermi-level crossing could be 
seen in ARPPES at overdoping where mobile bipolarons coexist with degenerate polarons 
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II. THE "FROHLICH-COULOMB" MODEL (FCM) 
A. Canonically transformed Hamiltonian 

Experimental facts tell us that any realistic description of high temperature supercon- 
ductivity should treat the long-range Coulomb and unscreened e-ph interactions on an equal 
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26] 



footing. In the past decade we have developed a "Frohlich-Coulomb" model (FCM) 
to deal with the strong long-range Coulomb and the strong long-range e-ph interactions in 
cuprates and other related compounds. The model Hamiltonian explicitly includes a long- 
range electron-phonon and the Coulomb interactions as well as the kinetic and deformation 
energies. The implicitly present large Hubbard U term prohibits double occupancy and 
removes the need to distinguish fermionic spins since the exchange interaction is negligible 
compared with the direct Coulomb and the electron-phonon interactions. 

Introducing fermionic, c n , and phononic, d ma , operators the Hamiltonian of the model is 
written as 



n^n' 



T(n - n')4c n / - 7^V c (n - - 



- u a g a (m - n)(e a ■ u m 

a,nm 

+ J^ 00 * (dlnadma + l/2) , (2) 
ma 

where T(n) is the hopping integral in a rigid lattice, e a is the polarization vector of the 
ath vibration coordinate, u m _ n = (m — n)/|m — n| is the unit vector in the direction from 
electron n to ion m, g a (m — n) is the dimensionless e-ph coupling function, and V c (n — n') 
is the inter-site Coulomb repulsion. g a (m — n) is proportional to the force / m ( n ) acting 
between the electron on site n and the ion on m. For simplicity, we assume that all the 
phonon modes are non-dispersive with the frequency u a , and include spin in the definition 
of n. We also use h = k B = c = 1. 

The phonon frequency dispersion can be readily included in the definition of all essential 
physical quantities such as the polaron level shift, mass and the polaron-polaron interaction 
using the quasi-momentum representation for phonons [3]. Generally it leads to a lighter po- 
laron compared with the nondispersive approximation. For example, comprehensive studies 
of the molecular Holstein model, in which the dispersive features of the phonon spectrum 
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are taken into account, found much lower values of the polaron mass compared with the 
non-dispersive model 271 ] . 

If the electron-phonon interaction is strong, i.e. the conventional e-ph coupling constant 
of the BCS theory is large, A > 1, then the weak-coupling BCS [28( and the intermediate- 



coupling Migdal-Eliashberg [29|, |30| approaches cannot be applied [3jJ. Nevertheless the 
Hamiltonian, Eq.(2), can be solved analytically using the "1/A" multi-polaron expansion 
technique jjj, if A = E p /zT(a) > 1. Here the polaron level shift is 

E p = ^2u a gl(n)(e a ■ u n ) 2 , (3) 

no: 

and zT(a) is about the half-bandwidth in a rigid lattice. The model shows a rich phase 
diagram depending on the ratio of the inter-site Coulomb repulsion V c and the polaron level 
shift E p [26]. The ground state of FCM is a polaronic Fermi liquid when the Coulomb 
repulsion is large, a bipolaronic high-temperature superconductor at intermediate Coulomb 
repulsions, and a charge-segregated insulator if the repulsion is weak. FCM predicts su- 
perlight polarons and bipolarons in cuprates with a remarkably high superconducting crit- 
ical temperature. Cuprate bipolarons are relatively light because they are inter-site rather 
than on-site pairs due to the strong on-site repulsion, and because mainly c-axis polarized 
optical phonons are responsible for the in-plane mass renormalization. The relatively small 
mass renormalization of polaronic and bipolaronic carries in FCM has been confirmed nu- 
merically using the exact Quantum Monte Carlo (QMC) {32], 3^], cluster diagonalization 



341 ] and variational [35] algorithms. 

The "1/A" expansion technique is based on the fact, known for a long time, that there is 
an analytical exact solution of a single polaron problem in the strong-coupling limit A — > 00. 
Following Lang and Firsov [36] we apply the canonical transformation e s diagonalising the 
Hamiltonian, Eq.(2). The diagonalisation is exact, if T(m) = (or A = 00). In the Wannier 
representation for electrons and phonons, 

S= g a (m - n)(e a • u m _ n )c| 1 c n (d| lia - d ma ). 

The transformed Hamiltonian is 

H = e~ s He s = - ^ + uo ^ \S mct d ma + - J + (4) 

n^n' ma ^ ' 



n^n' 



where for simplicity we take uo a = u. The last term describes the energy gained by polarons 
due to the e-ph interaction. The third term on the right-hand side is the polaron-polaron 
interaction, 

v(n - n') = V c (n - n') - V ph (n - n'), (5) 

where 



x 



V^(n-n') = 2u^2g a (xn-n)g a (m-ri) 

(Gq, ■ U m _ n )(GQ, • U m _ n '). 

The phonon-induced interaction V p h is due to displacements of common ions caused by two 
electrons. Finally, the transformed hopping operator <j nn / is given by 



Cnn' 



T(n — n') exp 



^ [g Q (m - n)(e Q ■ u r 



(6) 



g a (m - n')(e Q • u m 

— n')] i^lna <^mct)] • 



This term is a perturbation at large A. It accounts for the polaron and bipolaron tunnelling 
and high temperature superconductivity [3| . In particular crystal structures like perovskites, 
a bipolaron tunnneling could appear already in the first order in T(n) (see below), so that 
<r nn ' can be averaged over phonon vacuum, if the temperature is low enough, T uj. The 
result is 



t(n - n') = ({&nn')) ph = T(n - n') exp[-^ 2 (n - n')], 



(7) 



where 



x 



5f 2 (n-n') = ^2g a (m - n)(e a ■ u m _ n 

[g a (m - n)(e Q • u m _ n ) - g a (m - ri)(e a ■ u m _ n /)] . 

By comparing Eqs.(7) and Eqs.(3,5), the bandwidth renormalization exponent can be ex- 
pressed via E p and V p h as follows 



9 (n- n 



1 

CO 



E p -~V ph (n-n') 



(8) 



In zero order with respect to the hopping the Hamiltonian, Eq.(4) describes localised po- 
larons and independent phonons, which are vibrations of ions around new equilibrium po- 
sitions depending on the polaron occupation numbers. The phonon frequencies remain 



FIG. 2: Simplified model of the copper-oxygen perovskite layer 26]. 



unchanged in this limit. The middle of the electron band falls by the polaron level-shift E p 
due to a potential well created by lattice deformation. The finite hopping term leads to the 
polaron tunnelling because of degeneracy of the zero order Hamiltonian with respect to site 
positions of the polaron. 



B. Super light small bipolarons in FCM: root to room temperature superconduc- 
tivity 

Now let us consider in-plane bipolarons in a two-dimensional lattice of ideal octahedra 
that can be regarded as a simplified model of the copper-oxygen perovskite layer, Fig. 2 [26 ]. 
The lattice period is a = 1 and the distance between the apical sites and the central plane 
is h — a/2 = 0.5. For mathematical transperancy we assume that all in-plane atoms, both 
copper and oxygen, are static but apex oxygens are independent three-dimensional isotropic 
harmonic oscillators. 

Due to poor screening, the hole-apex interaction is purely coulombic, 

g a (m-n) = -, 

m — n\ z 

where a — x,y, z. To account for the experimental fact that ^-polarized phonons couple 
to the holes stronger than others we choose k x = K y = k z /^/2. The direct hole-hole 
repulsion is 

K(n-n') l ' 



v^ln-n'l 
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FIG. 3: Four degenerate in-plane bipolaron configurations A, B, C, and D . Some single-polaron 
hoppings are indicated by arrows [261 ] . 



so that the repulsion between two holes in the nearest neighbour (NN) configuration is V c . 
We also include the bare NN hopping T NN , the next nearest neighbour (NNN) hopping 
across copper T^nn and the NNN hopping between the pyramids T' NNN . 

The polaron shift is given by the lattice sum Eq.(3), which after summation over polar- 
izations yields 

E p = £ (r^n + T^Je) = 31.15«X (9) 

m \ i i / 

where the factor 2 accounts for two layers of apical sites. For reference, the Cartesian 

coordinates are n = (n x + l/2, n y + l/2, 0), m = (m x , m y , h), and integers. 

The polaron-polaron attraction is 

i .. 9 \ -v h 2 + (m — n') • (m — n) . 

V ph n - n' = 4w«* —T^ ,J 1 ,3 • 10 

/ — ' m — n'rlm — nr 

m 1 11 1 

Performing the lattice summations for the NN, NNN, and NNN' configurations one finds 
V p h = 1.23 E p , 0.80 E p , and 0.82 E p , respectively. As a result, we obtain a net inter-polaron 
interaction as vnn = V c — 1.23 E p , vnnn = % — 0.80 E p , v' NNN = ^| — 0.82^, and the 
mass renormalization exponents as g 2 NN = 0.38(E p /u), g 2 NNN = 0.60(Ep/uj) and {<]' nnn) 2 = 
0.59(E p /lu). 
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Let us now discuss different regimes of the model. At V c > 1.23 E p , no bipolarons are 
formed and the system is a polaronic Fermi liquid. Polarons tunnel in the square lattice 
with t = T NN exp(— 0.38E p /u) and t' = T NNN exp(—0.60E p /u) for NN and NNN hoppings, 
respectively. Since g 2 NNN ~ (q'nnn) 2 one can neglect the difference between NNN hoppings 
within and between the octahedra. A single polaron spectrum is therefore 

^(k) = -E p -2t'[cosk x + cos k y ] - Atcos(k x /2) cos(fe y /2). (11) 

The polaron mass is m* = l/(t + 2t'). Since in general t > t', the mass is mostly determined 
by the NN hopping amplitude t. 

If V c < 1.23 E p then intersite NN bipolarons form. The bipolarons tunnel in the plane 
via four resonating (degenerate) configurations A, B, C, and D, as shown in Fig. 3. In the 
first order of the renormalised hopping integral, one should retain only these lowest energy 
configurations and discard all the processes that involve configurations with higher energies. 
The result of such a projection is the bipolaron Hamiltonian 

H b = (V c - 3.23 E p ) + B \ B \ + C/Ci + D\d x ] (12) 

i 

-H Y^i A \Bi + Bid + Cf A + DlAi + H.c] 
i 

~t' Y\ A l*Bl + B ly C l + ClVxA + D LyAl + H.C.}, 

n 

where 1 numbers octahedra rather than individual sites, x = (1,0), and y = (0, 1). A Fourier 
transformation and diagonalization of a 4 x 4 matrix yields the bipolaron spectrum: 

E 2 (K) = V C - 3.23£ p ± 2t'[cos(K x /2) ± cos{K y /2)]. (13) 

There are four bipolaronic subbands combined in the band of the width 8t'. The effective 
mass of the lowest band is m** = 2/t'. The bipolaron binding energy is A w 1.23E P — V c . 
Inter-site bipolarons already move in the first order of the single polaron hopping. This 
remarkable property is entirely due to the strong on-site repulsion and long-range electron- 
phonon interactions that leads to a non-trivial connectivity of the lattice. This fact combines 
with a weak renormalization of t' yielding a superlight bipolaron with the mass m** oc 
exp(0.60 Ep/u). We recall that in the Holstein model m** oc exp(2E p /u)) [18|. Thus the 
mass of the Frohlich bipolaron in the perovskite layer scales approximately as a cubic root 
of that of the Holstein bipolaron. 
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At even stronger e-ph interaction, V c < 1.16-Ep, NNN bipolarons become stable. More 
importantly, holes can now form 3- and 4-particle clusters. The dominance of the potential 
energy over kinetic in the transformed Hamiltonian enables us to readily investigate these 
many-polaron cases. Three holes placed within one oxygen square have four degenerate 
states with the energy 2{V C — 1.23E P ) + V c j 'V2 — 0.80E P . The first-order polaron hopping 
processes mix the states resulting in a ground state linear combination with the energy 
E 3 = 2.71V C — 3.26-Ep — \/At 2 + 1' 2 . It is essential that between the squares such triads 
could move only in higher orders of polaron hopping. In the first order, they are immobile. 
A cluster of four holes has only one state within a square of oxygen atoms. Its energy is 
E 4 = 4{V C -1.23E P )+2(V C /V2-0.80E P ) = 5A1V C -6.52E P . This cluster, as well as all bigger 
ones, are also immobile in the first order of polaron hopping. We would like to stress that 
at distances much larger than the lattice constant the polaron-polaron interaction is always 
repulsive, and the formation of infinite clusters, stripes or strings is prohibited [371 ] . Hence 
the long-range Frohlich interaction combined with Coulomb repulsion might cause clustering 
of polarons into finite-size quasi-metallic mesoscopic textures. However analytical 37J] and 
QMC 38] studies of mesoscopic textures with lattice deformations and Coulomb repulsion 
show that pairs (i.e. bipolarons) dominate over phase separation since they effectively repel 
each other [3]. 

As shown above FCM is reduced to an extended Hubbard model with intersite attraction 
and suppressed double-occupancy in the limit of high phonon frequency uj > T(a) and large 
on-site Coulomb repulsion. Then the Hamiltonian can be projected onto the subspace of 
nearest neighbor intersite bipolarons, Fig. 3. In contrast with the crawler motion of on-site 
bipolaron, the intersite bipolaron tunnelling is a crab-like, so that its mass scales linearly with 
the polaron mass [m** « Am* on the staggered chain [26J). To study FCM for more realistic 
intermediate values of the electron-phonon coupling and phonon frequency the CTQMC 
algorithm 32j, ]39] has been recently extended to systems of two particles with strong electron- 
phonon interactions. We have solved the bipolaron problem on a stag gered ladder, triangular 
and anisotropic hexagonal lattices from weak to strong coupling 33J in a realistic parameter 
range where usual limiting approximations fail. The bipolaron to polaron mass ratio has 
been found about 2 in the weak coupling regime (A <C 1) as it should be for a large bipolaron 
4p| . In the strong-coupling, large phonon frequency limit the mass ratio approaches 4, in 
agreement with strong-coupling arguments given above. In a wide region of parameter space, 
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we find a bipolaron/polaron mass ratio of between 2 and 4 and a bipolaron radius similar 
to the lattice spacing. Thus the bipolaron is small and light at the same time. Taking 
into account additional intersite Coulomb repulsion V c does not change this conclusion. As 
V c increases the bipolaron mass decreases but the radius remains small, at about 2 lattice 
spacings. 

When bipolarons are small and pairs do not overlap, the pairs can form a Bose-Einstein 
condensate (BEC). Our CTQMC simulations show that with realistic values for the coupling 
constant, A ~ 1, and phonon frequencies, uo ~ T(a) one can avoid overlap of pairs and get 
the bose-condensation temperature T c about the room temperature. We believe that the 
following recipe is worth investigating to look for room-temperature superconductivity 33] : 
(a) The parent compound should be an ionic insulator with light ions to form high-frequency 
optical phonons, (b) The structure should be quasi two-dimensional to ensure poor screening 
of high-frequency c-axis polarized phonons, (c) A triangular lattice is required in combination 
with strong, on-site Coulomb repulsion to form the small superlight Crab bipolaron (d) 
Moderate carrier densities are required to keep the system of small bipolarons close to the 
dilute regime. 

PI 

There are strong arguments in favor of 3D bipolaronic BEC in cuprates [3J] drawn u sing 
parameter- free fitting of experimental T c with BEC T c 4JJ, unusual upper critical fields 
and the specific heat |43lj . and, more recently normal state diamagnetism 



Lorenz numbers 
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4611 . the normal state Nernst effect 
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44j, the Hall- 



481 ]. and the giant proximity 



effect (GPE) |49|] as discussed below. 



III. SOME NORMAL STATE PROPERTIES OF CUPRATES IN FCM 

A. Normal state Nernst effect and insulating-like in-plane resistivity 

The low-energy FCM electronic structure of cuprates is shown in Fig. 4 [50(. Polaronic 
p-holes are bound in lattice inter-site singlets (A) or in singlets and triplets (B) (if spins 
are included in Eq.(2)) at any temperature. Above T c a charged bipolaronic Bose liquid is 
non-degenerate and below T c phase coherence (ODLRO) of the preformed bosons sets in. 
The state above T c is perfectly " normal" in the sense that the off-diagonal order parameter 
(i.e. the Bogoliubov-Gor'kov anomalous average jF(r, r') = (ifji(r)ip^(r')) is zero above the 
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FIG. 4: Bipolaron picture of high temperature superconductors. A corresponds to a singlet oxygen 
intersite bipolaron, B is a triplet intersite bipolaron, A is the singlet bipolaron binding energy, J 
is the singlet-triplet exchange energy, and 2t is the bipolaron bandwidth 5C|. 



resistive transition temperature T c . Here ipi^(r) annihilates electrons with spin |, f at point 



In contrast with our bipolaronic (and with BCS) theory a significant fraction of research 
in the field of cuprate superconductors suggests a so-called phase fluctuation scenario, where 
jF(r, r') remains nonzero well above T c . I believe that the phase fluctuation scenario is im- 
possible to reconcile with the extremely sharp resistive transitions at T c in high-quality 
underdoped, optimally doped and overdoped cuprates. For example, the in-plane and out- 
of-plane resistivity of Bi — 2212, where the anomalous Nernst signal has been measured 
51] , is perfectly "normal" above T c , Fig. 5, showing only a few percent positive or negative 



magnetoresistance 



and out-of-plane [5 



52 1 , explained with bipolarons 53]. Both in-plane 54 



55 



56 



60 



53, 



58) 



611 ] resistive transitions of high-quality samples remain sharp in 
the magnetic field providing a reliable determination of the genuine H C 2(T). The preformed 
Cooper-pair (or phase fluctuation) model 62] is incompatible with a great number of ther- 



modynamic, magnetic, and kinetic measurements, which show that only holes (density x), 
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FIG. 5: In-plane (A) and out-of-plane (B) resistivity of 3 single crystals of Bi2Sr2CaCu20g 
showing no signature of phase fluctuations above the resistive transition. 



doped into a parent insulator are carriers both in the normal and the superconducting states 



of cuprates. The assumption 62J that the superfluid density x is small compared with the 
normal-state carrier density is also inconsistent with the theorem 6^ . which proves that 
the number of supercarriers at T = OK should be the same as the number of normal-state 
carriers in any clean superfluid. 

Recently we have described a number of unusual normal state properties in cuprates in a 
different manner as perfectly normal state phenomena. In particular, the bipolaron theory 
accounts for the anomalously large Nernst signal, the thermopower and the insulating-like 



in-plane low temperature resistance 



47 



48( as observed 



51 



64j, 



65, 



66|]. 



Thermomagnetic effects appear in conductors subjected to a longitudinal temperature 
gradient V X T in x direction and a perpendicular magnetic field in z direction. The transverse 
Nernst-Ettingshausen effect 63] (here the Nernst effect) is the appearance of a transverse 
electric field E y in the third direction. When bipolarons are formed in the strong-coupling 
regime, the chemical potential is negative. It is found in the impurity band just below 
the mobility edge at T > T c . Carriers, localised below the mobility edge contribute to 
the longitudinal transport together with the itinerant carriers in extended states above the 
mobility edge. Importantly the contribution of localised carriers of any statistics to the 
transverse transport is normally small [68j since a microscopic Hall voltage will only develop 
at junctions in the intersections of the percolation paths, and it is expected that these are 
few for the case of hopping conduction among disorder- localised states [69(. Even if this 
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contribution is not negligible, it adds to the contribution of itinerant carriers to produce 
a large Nernst signal, e y (T, B) = —E y /V X T, while it reduces the thermopower S and the 
Hall angle G. This unusual "symmetry breaking" is completely at variance with ordinary 



metals where the familiar "Sondheimer" cancelation [70| makes e y much smaller than Stan 
because of the electron-hole symmetry near the Fermi level. Such behaviour originates in 
the "sign" (or n p — n") anomaly of the Hall conductivity of localised carriers. The sign of 
their Hall effect is often opposite to that of the thermopower as observed in many amorphous 



semiconductors 



681 ] and described theoretically [71]. 



The Nernst signal is expressed in terms of the kinetic coefficients au and cn^ as 



_ &xx®-yx Gyx^-xx I -\ A\ 

e V = a 2 +a 2 ' W 

xx 1 xy 

where the current density is given by ji — a ijEj + otijV jT. When the chemical potential 
p is at the mobility edge, localised carriers contribute to the transport, so and can 
be expressed as af? 1 + a\, and afj* + a l ij, respectively. Since the Hall mobility of carriers 
localised below /z, a l yx , has the sign opposite to that of carries in the extended states above 
/i, cr^f, the sign of the off-diagonal Peltier conductivity a l yx should be the same as the sign 
of ot y x x . Then neglecting the magneto-orbital effects in the resistivity (since B -C 1 ) we 
obtain 



StanO = -^L_ „ pK f _ \a l xx \)(® ext ~ (15) 

and 



a 2 +a 2 

^ xx ' xy 



e y ^P« + Kl)- 5tan0 ' ( 16 ) 
where Q ext = /a xx , Q l = a l yx /a xx , and p = l/cr xx is the resistivity. 

Clearly the model, Eqs. (15,16) can account for a low value of StanB compared with a 
large value of e y in some underdoped cuprates 5l|, [65] due to the sign anomaly. Even in 
the case when localised bosons contribute little to the conductivity their contribution to 
the thermopower S = p(o> e xx — \a l xx \)) could almost cancel the opposite sign contribution 
of itinerant carriers [47|. Indeed the longitudinal conductivity of itinerant two-dimensional 
bosons, a ext oc J Q dEEdf(E)/dE diverges logarithmically when \x in the Bose-Einstein dis- 
tribution function f(E) = [exp((i? — /i)/T) — 1] _1 goes to zero and the relaxation time r is a 
constant. At the same time a e xx oc j dEEiE — p)df (E) / dE remains finite, and it could have 
the magnitude comparable with a l xx . Statistics of bipolarons gradually changes from Bose 
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FIG. 6: Normal state in-plane resistivity of underdoped Lai.94 Sro.o6Cu04 (triangles 



641 ]) as re- 



vealed in the field B = 12 Tesla and compared with the bipolaron theory, Eq.(19) (solid line). 

to Fermi statistics with lowering ener gy across the mobility edge because of the Coulomb 
repulsion of bosons in localised states [72] . Hence one can use the same expansion near the 
mobility edge as in ordinary amorphous semiconductors to obtain the familiar textbook re- 
sult S = SoT with a constant Sq at low temperatures [73(. The model becomes particularly 
simple, if we neglect the localised carrier contribution to p, and a xy , and take into account 



that o& oc B I p and Q ext oc B/p in the Boltzmann theory. Then Eqs. (15,16) yield 



and 



^tanB oc T/p 



e v (T,B)oc(l-T/T 1 )/p. 



(17) 



(18) 



According to our earlier suggestion [74| the insulating-like low-temperature dependence of 
p(T) in underdoped cuprates originates from the elastic scattering of nondegenerate itiner- 
ant carriers off charged impurities. We assume here that the carrier density is temperature 
independent at very low temperatures. The relaxation time of nondegenerate carriers de- 
pends on temperature as t oc T~ 1//2 for scattering off short-range deep potential wells, and 
as T 1 / 2 for very shallow wells jz^J]. Combining both scattering rates we obtain 



p = p [(T/T 2 ) 1 / 2 + (T 2 /T) 1 / 2 ]. 



(19) 



Eq.(19) with p = 0.236 mO-crn and T 2 = 44. 6K fits extremely well the experimental 
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insulating-like normal state resistivity of underdoped Lai.94 Sr .06CuO4 in the whole 



temperature range from 2K up to 50K, Fig. 6, as revealed in the field B = 12 Tesla 64 



ow- 



65] 



Another high quality fit can be obtained combining the Brooks-Herring formula for the 



3D scattering off screened charged impurities, as proposed in Ref.[75| for almost undoped 
LSCO, or the Coulomb scattering in 2D (r oc T) and a temperature independent scattering 
rate off neutral impurities with the carrier exchange [76J similar to the scattering of slow 
electrons by hydrogen atoms in three dimensions. Hence the scale T2, which determines the 
crossover toward an insulating behavior, depends on the relative strength of two scattering 
mechanisms. Importantly our expressions (17,18) for S*tan0 and e y do not depend on the 
particular scattering mechanism. Taking into account the excellent fit of Eq.(19) to the 
experiment, they can be parameterized as 

(TIT ) 3 / 2 

£tan9 = e K 1 , (20) 
1 + T/T 2 v ; 

and 

e y {T,B) = e — - , (21) 

where T\ and e are temperature independent. 

In spite of all simplifications, the model describes remarkably well both S"tanB and e y , 
Fig. 7, measured in Lai. 94 Sr .06CuO4 with a single fitting parameter, T\ = 50K using the 

17 



experimental p(T). The constant eo = 2.95 /xV/K scales the magnitudes of S"tanB and 
e y . The magnetic field B — 12 Tesla destroys the superconducting state of the low-doped 
La x 94 Sr .06CuO 4 down to 2K, Fig.6, so any residual superconducting order above 2K is 
clearly ruled out, while the Nernst signal, Fig. 7, is remarkably large. The coexistence of the 
large Nernst signal and a nonmetallic resistivity is in sharp disagreement with the vortex 
scenario, but in agreement with our model. Taking into account the field dependence of 
the conductivity of localised carriers, the pho non-drug effect, and their contribution to the 
transverse magnetotransport can well describe the magnetic field dependence of the Nernst 
signal {47I and improve the fit in Fig. 7 at the expense of the increasing number of fitting 
parameters. 



B. Hall-Lorenz number 



Recent measurements of the Righi-Leduc effect provides further evidence for real-space 
charged bosons preformed above T c 45|, |46[. The effect describes transverse heat flow re- 
sulting from a perpendicular temperature gradient in an external magnetic field, which is 
a thermal analog of the Hall effect. Using the effect the "Hall-Lorenz" electronic num- 
ber, L H = ie/hs) 2 K xy /(Ta xy ) has been directly measured 77] in YBa 2 Cu 3 OQ, 95 and 
YBa 2 Cu 3 6 . 6 since transverse thermal n xy and electrical a xy conductivities involve pre- 
sumably only electrons. The experimental Lh(T) showed a quasi-linear temperature de- 
pendence above the resistive T c , which strongly violates the WF law. Remarkably, the 
measured value of L H just above T c turned out precisely the same as predicted by the bipo- 
L = 0.15Lq, where Lq = tt 2 /3 is the conventional Sommerfeld value. The 



laron theory 



breakdown of the WF law revealed in the Righi-Leduc effect 



771 ] has been explained by a 



temperature-dependent contribution of thermally excited single polarons to the transverse 



magneto-transp ort 



45j. 



Surprisingly more recent measurements of the Hall-Lorenz number in single crystals of 
optimally doped YBa 2 Cu 2 6 ,g 5 and optimally doped and underdoped EuBa 2 Cu 3 O y led to 
an opposite conclusion [79]. The experimental Lu for these samples has turned out only 
weakly temperature dependent and exceeding the Sommerfeld value by more than 2 times 
in the whole temperature range from T c up to the room temperature. Following an earlier 



claim|80j Matusiak and Wolf |79| have argued that a possible reason for such significant 
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difference might be that Zhang et al. [77| used different samples, one for n xy and another 
one for a xy measurements, which makes their results for Lh inconsistent. 

Actually it has been shown 46j that there is no inconsistency in both Lh determinations. 
One order of magnitude difference in two independent direct measurements of the normal- 
state Hall-Lorenz number is consistently explained by the bipolaron theory |3j. The theory 
explains the huge difference in the Hall-Lorenz numbers by taking into account the difference 

n . n 

between the in-plane resistivity of detwinned [77[ and twinned [79J single crystals. It fits 
well the observed L#(T)s and explains a sharp Hall-number maximum [79J observed in the 
normal state of underdoped cuprates. 

In the presence of the electric field E, the temperature gradient VT and a weak magnetic 
field B|| z _L E and VT, the electrical currents in x,y directions are given by 

jx=a xx V x (u, - 2e0) + a xy Vy(u, - 2e</>) 
j y =a yy Vy(i2 - 2e0) + a yx V x (/j, - 2e<p) 

(22) 



and the thermal currents are: 



w x =c xx V x (/2 - 2e0) + c xy V y (/2 - 2e0) 

'^~d XX '^ X T ~\- d X y^ yT 

Wy=CyyVy(fl ~ 2e(j)) + Cy X V X (fi ~ 2C(j)) 
-\-dyy^ yT ~\~ (ly X ^ X T . 



(23) 



Here \x and are the chemical and electric potentials. 

Real phonons and (bi)polarons are well decoupled in the strong-coupling regime of the 
electron-phonon interaction [3| so the standard Boltzmann equation for the kinetics of renor- 
malised carriers is applied. If we make use of the r{E)~ approximation [81| the kinetic 
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coefficients of bipolarons are found as 45 ] 
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where 



(Q(E)) 



/ °° dEQ(E)ED b (E)df b / dE 



(24) 



/ °° dEED b (E)df b /dE ' 
D b (E) oc E d l 2 ~ l is the density of states of a <i-dimensional bipolaron spectrum, E = 
K 2 /(2m b ), g b = 2e/m b , and f b (E) is the equilibrium distribution function. Polaronic coef- 
ficients are obtained by replacing super/subscripts b for p, double elementary charge 2e for 
e and p for p/2 in all kinetic coefficients, and m b for 2m p in and Qj. The kinetic energy 
of bipolarons, E should be replaced by E + T*, where E = k 2 /(2m p ) is the polaron kinetic 
energy, and T* is half of the bipolaron binding energy (i.e. the pseudogap temperature in 
the theory |3j). 

The in-plane resistivity, p, the Hall number, Rh, and the Hall-Lorenz number, Lh are 
expressed in terms of the kinetic coefficients as p -1 = 2ea xx , Rh = a yx /2eB(a xx ) 2 , and 



^ \(dyxQ>x 



Cy X b xx )(lxx 



Cxx(bxx@>yx by X Qj XX y^ 



2J. (Xy X Qi xx 



(25) 



respectively, where a, b,c,d = aP + a b , IP + b b , c p + c b , d p + d b . 

The in-plane resistivity, the temperature-dependent paramagnetic susceptibility, and the 
Hall ratio have been described by the bipolaron model taking into account thermally ac- 



tivated single polarons 



82 



83 



84 



851 ] . The bipolaron model has also offered a simple 
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explanation of c-axis transport and the anisotropy of cuprates [24|, |60j, |6l|, |84J . The crucial 
point is that single polarons dominate in c-axis transport at finite temperatures because they 
are much lighter than bipolarons in c-direction. Bipolarons can propagate across the planes 
only via a simultaneous two-particle tunnelling, which is much less probable than a single 
polaron tunnelling. However, along the planes polarons and inter-site bipolarons propagate 
with comparable effective masses (section 2). Hence in the mixture of nondegenerate quasi- 
two-dimensional (2D) bosons and thermally excited 3D fermions, only fermions contribute to 



c -axis transport, if the temperature is not very low, which 



c -axis transport and to the huge anisotropy of cuprates 241 ] . 



eads to the thermally activated 



We have also shown [45J that by the necessary inclusion of thermally activated polarons, 
the model predicts a breakdown of the WF law with the small near-linear in temperature 
Hall-Lorenz number, as observed experimentally by Zhang et al. [ztJ (see Fig.8). Let us 
now show that the bipolaron model describes the contrasting observations of Ref. 3] as 
well, if the ratio of bipolaron and polaron mobilities, a = 2Tbm p /T p mb becomes relatively 
small. 

Both polaronic and bipolaronic carriers are not degenerate above T c , so the classical 
distribution functions, /& = yexp(— E/T) and f p = y 1 ^ 2 exp[— (E + T*)/T] are applied with 
y = exp(/i/T). The chemical potential is evaluated using 2nb + n p = x/vo, where x is the 
number of itinerant holes in the unit cell volume vq not localised by disorder. The bipolaron 
density remains large compared with the polaron density in a wide temperature range, so 
that ritVo ~ x/2 and y ~ 7tx/(mba 2 T) for quasi-2D bipolarons. Then the atomic density of 
3D polarons is u p vq = Tm p a 2 exp(— T* /T)(xm p /2ir 2 mb) lj/2 (a is the lattice constant). The 
ratio (3 = n p /2nb remains small at any pseudogap temperature T* and any relevant doping 
level x > 0.05, (3 ~ Texp(— T*/T)(18m p /7r 2 a;mft) 1 / 2 /jy <C 1, if the temperature T is small 
compared with the polaron bandwidth W = 6/m p a 2 . Hence, if the mobility ratio a is of the 
order of unity, both longitudinal and transverse in-plane magnetotransport is dominated by 
bipolarons, which explains a remarkably low Lh in high-quality detwinned crystals used in 
Ref.Q, Fig.8. 



On the other hand, twinned crystals used in Ref.[79| had the in-plane resistivity several 
times larger than those of Ref. 77| presumably resulting from twin boundaries and long 
term aging. The twin boundaries and other defects are strong scatterers for slow 2D bipo- 
larons, while lighter quasi-3D polarons are mainly scattered by real optical phonons, which 
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FIG. 8: The Hall-Lorenz number Lh in underdoped twinned EuBa2Cu^O^^ (circles) 



79f ] com- 



pared with the theory, Eq.(28) when a <C 1 (upper line), and the significantly different Hall-Lorenz 
number in detwinned YBa2Cu 3 OQ,^ (triangles) \l\ described by the same theory [2| with a mod- 
erate value of a = 0.44 (lower line). 



are similar in all crystals. Hence one can expect that a becomes small in twinned crystals 
of Ref. {to!]. If the condition a 2 <C (3 is met, then only polarons contribute to the transverse 
electric and thermal magnetotransport. It explains about the same thermal Hall conduc- 
tivities (n xy 2.5 x 10~ 3 .B W/Km at T=100K) dominated by polarons in both crystals of 
YBa 2 Cu 3 O e ,g 5 used in Ref. 77j and in Ref. [79(, and at the same time a substantial differ- 
ence of their electrical Hall conductivities, a xy , as bipolarons virtually do not contribute to 
a xy in the twinned samples. 

To arrive at simple analytical results and illustrate their quantitative agreement with 
the experiment jjS] let us assume that a 2 <C j3, but a > (3, and neglect an energy depen- 
dence of the transport relaxation rates of all carriers. In such conditions bipolarons do not 
contribute to transverse heat and electric flows, but determine the in-plane conductivity. 
Kinetic responses are grossly simplified as 



P = 2e^ b (26) 

a* = ^ = ^, (27) 

exa z mi 

L H = 4.75 + 3T*/T + (T*/T) 2 . (28) 
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As in the case of a 2 > /3, discussed in Ref. 45[], the recombination of a pair of polarons 
into bipolaronic bound states at the cold end of the sample results in the breakdown of 
the WF law, as described by two temperature-dependent terms in Eq.(28). The breakdown 
is reminiscent of the one in conventional semiconductors caused by the recombination of 
electron- hole pairs at the cold end 8J|. However, the temperature dependence and the value 
of Lh{T) turn out remarkably different. When a 2 -C /3, The Hall-Lorenz number is more 
than by an order of magnitude larger than in the opposite regime, a 2 > f3. It increases with 
temperature lowering rather than decreases fitting well the experimental observation [3] in 
twinned underdoped single crystals of EuBa^CuzO^,^ with T* = 100K, Fig.8. Hence by 
varying the bipolaron to polaron mobility ratio, a, the model accounts for qualitatively 
different behaviours of Lh(T) in twinned and detwinned cuprates. 



C. Normal state diamagnetism 



A number of experiments (see, for example, 66 



86 



87 



89 



90] and references therein) , 



including torque magnetometries, showed enhanced diamagnetism above T c , which has been 
explained as the fluctuation diamagnetism in quasi-2D superconducting cuprates (see, for 
example Ref. 88(). The data taken at relatively low magnetic fields (typically below 5 
Tesla) revealed a crossing point in the magnetization M(T, B) of most anisotropic cuprates 
(e.g. Bi - 2212), or in M(T,B)/B^ 2 of less anisotropic YBCO [STJ. The dependence of 
magnetization (or M/B 1 ^ 2 ) on the magnetic field has been shown to vanish at some char- 
acteristic temperature below T c . However the data taken in high magnetic fields (up to 30 
Tesla) have shown that the crossing point, anticipated for low-dimensional superconductors 
and associated with superconducting fluctuations, does not explicitly exist in magnetic fields 
above 5 Tesla 



SSI- 



Most surprisingly the torque magnetometery 86j, [89J uncovered a diamagnetic signal 
somewhat above T c which increases in magnitude with applied magnetic field. It has been 
linked with the Nernst signal and mobile vortexes in the normal state of cuprates 66]. 
However, apart from the inconsistences mentioned above, the vortex scenario of the normal- 
state diamagnetism is internally inconsistent. Accepting the vortex scenario and fitting 
;he magnetization data in Bi — 2212 with the conventional logarithmic field dependence 
66| . one obtains surprisingly high upper critical fields H C 2 > 120 Tesla and a very large 
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Ginzburg-Landau parameter, k — A/£ > 450 even at temperatures close to T c . The in- 
plane low-temperature ma gnet ic field penetration depth is A = 200 nm in optimally doped 
Bi — 2212 (see, for example 9JJ). Hence the zero temperature coherence length £ turns out to 
be about the lattice constant, £ = 0.45nm, or even smaller. Such a small coherence length 
rules out the "preformed Cooper pairs" 62], since the pairs are virtually not overlapped 
at any size of the Fermi surface in Bi — 2212 . Moreover the magnetic field dependence of 
M(T, B) at and above T c is entirely inconsistent with what one expects from a vortex liquid. 
While —M(B) decreases logarithmically at temperatures well below T c , the experimental 



curves 



66 



86 



89| clearly show that —M(B) increases with the field at and above T c , just 



opposite to what one could expect in the vortex liquid. This significant departure from the 
London liquid behavior clearly indicates that the vortex liquid does not appear above the 
resistive phase transition [861 ] . 



Some time ago we explained the anomalous diamagnetism in cuprates as the Landau 
normal-state diamagnetism of preformed bosons 92|. The same model predicted the unusual 



upper critical field 42| observed in many superconducting cuprates 
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More recently the model has been extended to high magnetic fields taking into account the 



44] 



magnetic pair-breaking of singlet bipolarons and the anisotropy of the energy spectrum 
When the magnetic field is applied perpendicular to the copper-oxygen plains the quasi-2D 
bipolaron energy spectrum is quantized as E a = oj{n + 1/2) + 2t c [l — cos(K z d)), where a 
comprises n — 0, 1, 2, ... and in-plane K x and out-of-plane K z center-of-mass quasi- moment a, 
uj = 2eB / ^Jm**m** , t c and d are the hopping integral and the lattice period perpendicular to 
the planes. We assume here that the spectrum consists of two degenerate branches, so-called 
"x" and "y" bipolarons as in the case of apex intersite pairs [25] with anisotropic in-plane 
bipolaron masses m** = m and m** ~ 4m. Expanding the Bose-Einstein distribution 
function in powers of exp[(/i — E)/T] with the negative chemical potential \i one can after 
summation over n readily obtain the boson density 

ixd ^— j' 1 - exp(-ujr/l ) 

and the magnetization, 

2eT A , f2t r r 



(29) 



M(T,B) = -n b , b+ 2 ^J2l ^Jf) 

r=l ^ ' 



X 



(30) 



exp[(^ - u/2 - 2t c )r/T] (\ 
1 — exp(— Lur/T) \r 



ti;exp(— ujr/T) 
ksTW — exp(— ur/T)] 
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FIG. 9: Diamagnetism of optimally doped Bi-2212 (symbols) [66] compared with magnetization of 
CBG {44] near and above T c (lines). 



Here //& = e/ ^/m**m** and Iq(x) is the modified Bessel function. At low temperatures 
T — > Schafroth's result 9j| is recovered, M(0, 5) = —n b jj, b . The magnetization of charged 
bosons is field-independent at low temperatures. At high temperatures, T ^> T c the chemical 
potential has a large magnitude, and we can keep only the terms with r = 1 in Eqs. (85,86) 
to obtain M(T,B) = — nb/if,u;/(6T) at T ^> T c ^> u, which is the familiar Landau or- 
bital diamagnetism of nondegenerate carriers. Here T c is the Bose-Einstein condensation 
temperature T c = 3.31(n 6 /2) 2 / 3 / {m**m**m**) l / z , with m c = l/2\t c \d 2 . 

Comparing with experimental data one has to take into account a temperature and 
field depletion of singlets due to their thermal excitations into spin-split triplets and single 
polaron states, Fig.4. If J < A/2, triplets mainly contribute to temperature and field 
dependencies of the singlet bipolaron density near T c , n&(T, B) = n c [l — ar — (B/B*) 2 ]. 
Here a = 3(2n c t)~ 1 [J(e J / Tc - l)" 1 - T c ln(l - e" J / Tc )], ft B B* = (2T c n c t) l l 2 sinh( J/2T C ), 
fis ~ 0.93 x 1CT 23 Am 2 is the Bohr magneton, n c is the density of singlets at T = T c in zero 
field, t = T/T c — 1, and 2t is the triplet bandwidth. A triplet contribution to diamagnetism 
remains negligible compared with the singlet diamagnetism if |1— T/T c \ J/T c . Also bosons 
localised in a random potential contribute to the diamagnetism. However their contribution 
remains small compared with the extended carrier diamagnetism, if the localization energy 
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is large compared with T. As a result, Eq.(30) fits remarkably well the experimental curves 
in the critical region of optimally doped Bi-2212, Fig.9, with n c fit> = 2100A/m, T c = 90K, 
a = 0.62 and B* = 56 Tesla, which corresponds to the singlet-triplet exchange energy 
J ss 20K. 



D. Giant proximity effect 

Several groups reported that in the Josephson cuprate SNS junctions supercurrent can 
run through normal iV-barriers thicker than 100 nm in a strong conflict with the standard 
theoretical picture, if the barrier is made from non-superconducting cuprates. Using an 



advanced molecular beam epitaxy, Bozovic et al. 94j proved that this giant proximity effect 
(GPE) is intrinsic, rather than extrinsic caused by any inhomogeneity of the barrier. Hence 
GPE defies the conventional explanation, which predicts that the critical current should 
exponentially decay with the characteristic length of about the coherence length, which is 
£ < 1 nm in the cuprates. 

Here I show that the effect can be broadly understood as the Bose-Einstein condensate 
(BEC) tunnelling into a cuprate semiconductor. As mentioned in section 1, the chemi- 
cal potential /i remains in the charge-transfer gap of doped cuprates like L^-zSr^CuO^ 
j?! because of bipolaron formation. The condensate wave function, ip(Z), is described by 
the Gross-Pitaevskii (GP) equation. In the superconducting region, Z < 0, near the SN 
boundary, Fig. 10, the equation is 

J-^ = [Vm)? - mZ) , (31) 

where V is a short-range repulsion of bosons, and m** is the boson mass along Z. Deep 
inside the superconductor ^(Z)! 2 = n s and \i = Vn s , where the condensate density n s is 
about x/2, if the temperature is well below T c of the superconducting electrode (the in-plane 
lattice constant a and the unit cell volume are taken as unity). 

The normal barrier at Z > is an underdoped cuprate semiconductor above its transition 
temperature, where the chemical potential \x lies below the bosonic band by some energy e, 
Fig. 10. One obtains for quasi-two dimensional bosons 3] 

e = -Tln(l -e~ To/T ), (32) 
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FIG. 10: BEC order parameter at the SN boundary for fi = 1.0,0.1,0.01 and ^ 0.001 (upper 



curve . 



where T = irx'/m**, m** is the in-plane boson mass, and x' < x is the doping level of the 
barrier. Then the GP equation in the barrier can be written as 

l d 2 4j{z) 



2m** dZ 2 



[VmZ)\* + eW(Z). (33) 



Introducing the bulk coherence length, £ = 1/ 1 (2m**n s V) l l 2 and dimensionless f(z) 
i)(Z) I 'nj 2 \ fx = e/n s V, and z = Z/£ one obtains for a real f(z) 



dz 2 

if z < 0, and 

d 2 f 



f ~ f, (34) 



r+fif, (35) 



dz 2 

if z > 0. These equations can be readily solved using first integrals of motion respecting the 
boundary conditions, /(— oo) = 1, and /(oo) = 0, 

f z = -(1/2 + f/2 - / 2 ) 1/2 , (36) 



and 

- = -(/V 2 + /72) 1/2 , (->7) 



dz 
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for z < and z > 0, respectively. The solution in the superconducting electrode is given by 



f(z) = tanh 



1/2 2V2(1 + £)V2 + 1 

2 z + U.51n 2l/2(i + ~ )1/2 _ 1 j 



(38) 



It decays in the close vicinity of the barrier from 1 to f(0) = [2(1 + //)] 1//2 in the interval 
about the coherence length £. On the other side of the boundary, z > 0, it is given by 

= sinh^/I 1 / 2 + ln[2(/i(l + /i)) 1 / 2 + (1 + 4/i(l + /i)) 1/2 ]}' ^ 
Its profile is shown in Fig. 10. Remarkably, the order parameter penetrates into the normal 
layer up to the length Z* « (/i) -1 / 2 ^, which could be larger than £ by many orders of 
magnitude, if fl is small. It is indeed the case, if the barrier layer is sufficiently doped. 
For example, taking x' = 0.1, c-axis m** = 2000m e , in-plane m** = 10m e jsj], a = 0.4 
nm, and f = 0.6 nm, yields T w 140 K and (/i)" 1/2 w 5000 at T = 10K. Hence the 
order parameter could penetrate into the normal cuprate semiconductor up to more than 



a thousand coherence lengths as observed 94J. If the thickness of the barrier L is small 



compared with Z*, and (/x) 1 ^ 2 1, the order parameter decays following the power law, 
rather than exponentially, 

m - (40) 



Hence, for L < Z*, the critical current should also decay following the power law 49(. On 



the other hand, for an undoped barrier \x becomes larger than unity, jl oc ln(m**T '/nx') — > oo 



for any finite temperature T when x' — > 0, and the current should exponential 
the characteristic length smaller that £, as is experimentally observed as well 

IV. DISCUSSION 



y decay with 

3- 



A possibility of real-space pairing, as opposed to the Cooper pairing, has been the subject 
of many discussions, particularly heated over the last 20 years after the discovery of high 
temperature superconductivity in cuprates. The first proposal for high temperature super- 
conductivity, made by Ogg Jr in 1946 [95J], already involved real-space pairing of individual 
electrons into bosonic molecules with zero total spin. This idea was further developed as 
a natural explanation of conventional superconductivity by Schafroth and Butler and Blatt 



961 ] . However, with one or two exceptions, the Ogg-Schafroth picture was condemned and 



practically forgotten because it neither accounted quantitatively for the critical behavior of 
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conventional superconductors, nor did it explain the microscopic nature of attractive forces 
which could overcome the Coulomb repulsion between two electrons constituting a pair. The 

individual electron pairs became fully transparent when 



failure of the 'bosonic' picture o 
Bardeen, Cooper and Schrieffer 28j] proposed that two electrons in a superconductor were 
indeed correlated, but at a very large distance of about 10 3 times of the average inter-electron 
spacing. 

Highly successful for low-Tc metals and alloys the BCS theory has led many researchers to 
believe that novel high-temperature superconductors should also be "BCS-like". However, 
the Ogg-Schafroth and the BCS descriptions are actually two opposite extremes of the same 
electron-phonon (e-ph) interaction. Indeed by extending the BCS theory towards the strong 
interaction between electrons and ion vibrations, a charged Bose gas (CBG) of tightly bound 



small bipolarons was predicted by us 



181 ] with a further prediction that high T c should exist 



in the crossover region of the e-ph interaction strength from the BCS-like to bipolaronic 



superconductivity 



19]. 



However, for very strong electron-phonon coupling, polarons become self-trapped on a 
single lattice site. The energy of the resulting small polaron is given as E p = —XzT(a), 
where A is the electron-phonon coupling constant, T(a) is the hopping parameter and z 
is the coordination number. Expanding about the atomic limit in small T(a) (which is 
small compared to E p in the small polaron regime, A > 1) the polaron mass is computed 
as m* = mo exp( / -fz\/hu>) , where u is the frequency of Einstein phonons, mo is the rigid 
band mass on a cubic lattice, and 7 is a numerical constant. For the Holstein model, which 
is purely site local, 7 = 1. Bipolarons are on-site singlets in the Holstein model and their 
mass m*^ appears only in the second order of T(a) 18j scaling as m*^ oc (m*) 2 in the limit 
ftijj ^> A , and as m# oc (m*) 4 in a more realistic regime hv <C A [97|]. Here A = 2E P — U is 
the bipolaron binding energy, and U is the on-site (Hubbard) repulsion. Since the Hubbard 
U is about 1 eV or larger in strongly correlated materials, the electron-phonon coupling 
must be large to stabilize on-site bipolarons and the Holstein bipolaron mass appears very 
large, m^/mo > 1000, for realistic values of phonon frequency. 

This estimate led some authors to the conclusion that the formation of itinerant small 
polarons and bipolarons in real materials is unlikely 98] , and high-temperature bipolaronic 



superconductivity is impossible 



991 ] . However, one should note that the Holstein model is 



an extreme polaron model, and typically yields the highest possible value of the (bi)polaron 
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mass in the strong coupling limit. Many advanced materials with low density of free carriers 
and poor mobility (at least in one direction) are characterized by poor screening of high- 
frequency optical phonons and are more appropriately described by the long-range Frohlich 



electron-phonon interaction 



251 ] . For this interaction the parameter 7 is less than 1 (7 « 0.3 



on the square lattice and 7 ~ 0.2 on the triangular lattice), reflecting the fact that in a hop- 
ping event the lattice deformation is partially pre-existent. Hence the unscreened Frohlich 
electron-phonon interaction provides relatively light small polarons, which are several or- 
ders of magnitude lighter than small Holstein polarons, which is now confirmed by several 
numerical studies as discussed in section 2. 

This unscreened Frohlich interaction combined with on-site repulsive correlations can bind 
holes into superlight intersite bipolarons (section 2). Experimental evidence for exceptiona 
strong e-ph interactions is now so overwhelming that the bipolaronic charged bose gas^ 
could be a feasible alternative to the BCS-like scenarios of cuprates. While some authors 
have dismissed any real-space pairing, advocating a collective pairing into incoherent Cooper 
pairs in the momentum space at some high temperature T* > T c , I argue that the most 
likely scenario is a true 3D Bose-Einstein condensation at T c of real-space bipolarons. Our 
bipolaron theory predicted such key features of cuprate superconductors as anomalous upper 
critical fields, spin and charge pseudogaps, and anomalous isotope effects later discovered 
experimentally. The theory explained normal state kinetics, high T c values, and specific heat 
anomalies of cuprates (see [3] and references therein). 

Here I have reviewed normal-state diamagnetism, the Nernst, thermal Hall and giant 
proximity effects as strong evidence for real-space pairing and 3D BEC in cuprates. 
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P.E. Kornilovitch, and J.H. Samson for illuminating discussions. The work was supported 
by EPSRC (UK) (grant no. EP/C518365/1). 
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